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1. Born's Reciprocity Principle 

One cannot help but be struck by the way that numerous structures in physics look the 
same under the simultaneous substitution between momentum p and position r, 

p^fer, r ^ -p/6, (1) 

where 6 is a scale with the dimensions of M/T. This applies to the classical Poisson 
brackets {rj,pj} = 5ij, the quantum commutator brackets = ifiSij, and the form 

of the Hamiltonian equations (classical or quantum), rj = dH/dpi, pi = —dH/dvi and 
of the angular momenta, Lij = riPj — VjPi. It leads to the concept of phase-space, Fourier 
transforms and uncertainty relations. The conjugacy between space and momentum is 
extensible to energy and time, this being required for special relativity. However these 
observations do not presume that Hamiltonians are invariant under transformation (1). 

Born's reciprocity principle [H [2] goes one stage further and assumes that all 
physical equations of motion are invariant and not just covariant under such conjugacy 
transformations, so that 

H{r,p)=H{-p/b,br). (2) 

At first sight this seems a patently absurd idea for anything but oscillators and even 
there it seems quite silly because it leads to fixed frequency for all vibrations, determined 
by the value of b. For these reasons and for its failure in accounting for the observed 
particle masses [Ij the principle has naturally fallen into disrepute and has never been 
taken seriously by physicists. There are also some fundamental philosophical objections 
to the idea, which will be mentioned later. In spite of these very valid criticisms, 
we wish to explore the principle and see if we can determine a non-zero value of b 
(which is probably tiny indeed and tied to cosmic scales). The incentive/reason why 
we wish to entertain the chance that eq.(2) may be valid arises also from vibrations. 
The point is that any oscillator contains some measure of anharmonicity; how this is 
manifested depends on the physical context, but we can be sure that the linear force and 
its associated potential V cannot increase without end. For instance we might suppose 
that in reality, for mass m, the true potential is V{r) ~ (mco'^r^/2) exp(— cr^), where c 
is a small parameter which sets the distance at which anharmonicity kicks in, and u is 
the natural rotational frequency for displacements that are not excessive. If we attempt 
to make H reciprocity invariant, we may arrive at 

2H = {p^ + b^r^)/m + mcu^p exp(-cr^) + (p^/b^) exp(-cp^/6^)]. 

It follows that if b is miniscule on ordinary momentum scales, the dangerous last term 
is minute, as is the correction to the kinetic energy, so the standard picture prevails for 
small or moderate r. This example indicates we have no right to be so dismissive of 
Born's principle. 

If we view the reciprocity substitution as the transformation 

( 6r ) ^ ( -1 ) ( 6r ) ' 
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we can think of another reciprocity transformation which also leaves most of our physical 
structures intact, namely 




ra im. 



In this paper we will consider a non-relativistic potential which is fully established 
on the large-scale (Newtonian gravity) and on the small scale (atomic physics), namely 
1/r; it has its roots in graviton and photon exchange and has been thoroughly 
studied over the centuries! We shall explore how Born's principle affects it. The 
immediate question is how to make V{r) oc 1/r compatible with (1). The substitution, 
1/r — i> 1/y^r^ + /h^ can be rejected outright as it would enormously enhance the 
velocity dependence for small 5, in fact ridiculously so. However a more reasonable 
alternative is 1/r (l/r + 6/p), since the last term fades out as 6 — > 0. There 
may exist other choices for making V compatible with Born's principle, but they are 
probably less natural than the proposal: 

where a signifies the interaction strength {Ze"^ /Aneo for a hydrogenic atom or GMm 
for gravity). At this stage we will take the sign of h positive even though it hails from 
1 / v^i? and we are not entirely sure about the sign of the root because the underlying 
dynamics (and associated field) is unclear. The choice of sign will become firmer in 
section 4 but it is fully consonant with the second form of reciprocity substitution. In 
short, for what it's worth, (3) will be our object of study It has the curious, if not 
dubious, feature that for any finite energy and r 7^ separation the speed can never 
vanish; the minimum v may be very small, being determined by the tiny constant 6, so 
it connotes particle restlessness (like zitterbewegung) even at the classical level! The 
only possibility for the speed to vanish is when the singularity at r = is reached. Since 
r stands for the relative distance between the test body and the centre of influence, it 
makes more sense to reinterpret h = mQ where m is the reduced mass and regard Q 
as the truly universal constant. Doing so will ensure that the modified kinetic energy 
assumes the same form no matter how coordinates are chosen: 

pj/mi + pl/m2 + mi^Vi + m2^'^rl = P^/M + p^/m + MVL^R^ + mVL^r'^. 

Here R is the centre of mass location, M = mi + m,2 is the total mass and P is the 
total momentum. Because of P, R dependence this does not mean that H is translation 
invariant any more than the usual Coulomb Hamiltonian is invariant under boosts; it 
is only covariant under those transformations. Therefore this represents a philosophical 
problem for relativity (even Galilean) and we shall worry about it later. 

I If a Yukawa potential exp(— /ir)/r is modified to obey reciprocity, the additional term is 
6exp(— and it becomes negiglible for small b; thus nuclear physics is unlikely to be affected. 
§ Such a Hamiltonian in turn spawns a strange-looking Lagrangian when it is expressed in terms of 
position and velocity. 
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The paper is set out as follows. Section 2 discusses the classical problem and 
trajectories as 6 is varied. Not unexpectedly we find that the distorted orbits precess 
around the force centre and we determine the rate for small 6; this is much like general 
relativistic corrections to Newtonian gravity. Section 3 deals with the quantum version; 
we find the change in energy levels to first order in h by the variational method and 
perturbation theory — which happen to agree with one another. In this way we set limits 
on the value of 6 so as not to disturb experimental atomic results, namely h < 10~^^ 
kg/s, a rather weak conclusion. More stringent limits come by looking at galactic 
scales, where the 1/v term can influence rotation rates profoundly. Section 4 contains 
our investigations of (3) for galaxies (possessing a supermassive black hole at the centre); 
there we find that the 1/v term can yield a velocity which at first increases linearly from 
the centre and then steadies out to a constant value, before rising again due to the effect 
of the harmonic accelaration. Our conclusions and continued worries with Born's 
principle end the paper in Section 5. 



2. Classical motion 



The equations of motion arising from the non-relativistic expression (3) are 

r = p(l/m + 6q;/p^), p = -r(6^/m + a/r^), (4) 

with H — E, the conserved "energy". This means there is a cubic relation between 
momentum and speed, force and displacement. For positive h the speed/momentum 
can never vanish, and if a > (an attractive interaction) the force cannot disappear 
either. For the purpose of the ensuing analysis we shall take 6 to be a small and positive 
quantity; many of the results can be continued to negative h without endangering the 
steps and we shall anyway be expressing the (small) precessions of the orbit to first 
order in h where the sign is somewhat irrelevant. 

It is quite difficult to solve these equations in general but it helps to use rotational 



invariance of (4) and the conserved angular momentum 
Hamiltonian in the usual way: 

02 



|r X p| to simplify the 



H 



1 

2m 



p. 



— a 



E: 



Pr = P • r 



(5) 



/p2 + £2/^2 

(Because of reciprocity we could equally well have used = r ■ p in place of Pr and p 
instead of br. However the form (5) is more familiar and this is the framework we shall 
adopt.) In consequence the radial equations are 



Pr_ 

m 



1 + 



abm 



(p2+£2/r2)3/2 



Pr 



m 



a 



+ 



[r'^Pr 



£2^3/2" 



(6) 



For determining the trajectory in the orbital plane, remember that the rate of change 
of azimuthal angle is ip = i/mr^, so the trajectory equation is obtained by integrating 

dip _ I 
dr r'^Pr[l + (xbm/ {i"^ /r"^ + p'^Y/'^y 



(7) 
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E (r,pr=0) 
20 




Figure 1. Dependence of energy on radius when the radial velocity vanishes. 



in which has to be ehminated in terms of E via (5). This is a hard problem for 
general h so we shall first turn to Mathematica for some elucidation of the motion and 
orbits. 

By looking at the E'-contours in phase space (r — pr) one finds that for £ 7^ the 
trajectories at lower energy are bounded with = r = at perigee or apogee. Absolute 
minima of E arise when dE/dr = or where 



h'^r a ah 
m r'^ i 



0. 



This cubic equation in r is readily solved and it has three real roots in the region of 
interest, 



am 



± 



To 



with corresponding extremal energy values 



'21 
am 



+ 2eb], Eo^ M 




(8) 



(9) 



2m V -^^ 

£■0 is an unstable maximum, while E± are equal value stable minima, and this is best 
revealed by plotting E{r, Pr — 0) against r in Figure 1. The phase space portrait is drawn 
in Figure 2 and the corresponding trajectory is depicted in Figure 3 over a timespan of 
two seconds. In all diagrams we have assumed unit mass and taken exaggerated values, 
£ — l,a — 10,b = 1 in order to emphasize the misshapen and precessing orbit. 

It is apparent that for 6 7^ orbits become distorted (sometimes very pronouncedly) 
from Kcplcrian ellipses and precession occurs. It is of interest to work out the 
precessional rate for model (3) and small b when the distortions/changes are tiny too. 
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Figure 3. Trajectory when E = —21. 
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To do so wc make the standard change of variable u — 1/r and expand (5) and (3) to 
first order in b. Thus 

-1 



dip__i_ 

du Pr 



1 + 



abm 



+ p2)3/2 



Pr 



1 - 



abm 



(£%2 + p2)3/2 



(10) 



and 



p — 2m{E + Qii)p — 2mnb + b u =0; p = Jp^ + i^u"^ 



The root we need is Pr ^ ^j2m{E + au) — Pu^, so expanding around that the trajectory 
equation (10) simphfies to 



dip 
du 



^2m{E + au) - i^u^ 
1 



1 - 



abm 



\J{U- Ui){u2 - U) 



1 



[2m{E + au)]^/\ 
abmi^ 

[(Ui + U2)u - UiU2f/^ 



(11) 



where Ui = l/ri,0i = at apogee and U2 = l/r2,(j)2 = •/r/2 at perigee - the turning 
points in r. The first term on the right of (11) produces the usual elliptical orbit answer 

2u — {ui + U2) + {ui — U2) cos{2if) 

and the integration of the second term is connected with the precession. Over an orbit 
the additional change in azimuth is 

- Aipb = abmf / , 77 ■ : tttt^ 

■^"1 ^J{U - Ul) {U2 - U) ii'^l + '^2)U - UiU2]^/^ 
; r del) 

Jo Uiul 



— abmi 



+ ul) + {uj - ul) COS0])3/2 



'^^-E{1 - ul/ul) = 2abmerir2E{l - rj/r^). 



(12) 



Here E{k) is the complete elliptic integral, which for small argument behaves as 
E(k) = {7r/2){l — k/4: — . . .) as k — > 0. Hence for small eccentricity (ri ~ r2 ~ a), 
the change in azimuth is 

Aipb ~ -27rabme^a^, (13) 

to a good approximation, where a is the semimajor axis. In the following section we shall 
set limits on b (for electrons at least) such that accurate atomic physics experiments are 
not substantially disturbed. 



3. Quantum mechanical considerations 

We will now be dealing with the operator version of (3) for hydrogenic atoms when 
a — Ze^/ATreQ. and 



H — — a 

2m 



+ 



(14) 
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It is evident that we are dealing with a tricky problem due to the last term, connected 
with 1/|-P|, even for eigenf unctions of angular momentum Ye^^{9^ (p) so £{£ + l)h'^ 

in the Schrodinger equation. 

Any serious attempt to try to solve (14) needs an interpretation of when i = 0. 
As PrR{r) —>■ {—ih/r)d{rR{r))/dr, when acting on the radial part R{r) of the wave 
function, a reasonable definition of the inverse is the indefinite integral: 

■^Rir) --- / dr'r'Rir') (15) 
nr Jr 

for solutions where rR{r) — > as r — > oo. One readily checks that Pj.P'^'^ — 
P~^Prip — ip for normalizable wave functions. It is even true that P~^ . exp{ikr) / r — 
(l/hk). exp{ikr)/r for outgoing waves, giving extra credence to the interpretation (15). 

In the event we have not succeeded in obtaining a complete solution of u{r) — ri/j{r) 
of type P(r) exp(— «;r — hr"^ /2%) to the equation for the radial wavefunction 

fi? d^u d (h'^r'^u au\ iabu du 
+ ^\^ -^ = E—, 



2m dr^ dr \ 2m r J h dr^ 

because of the troublesome last term on the lefthand side. No matter. One can still 
obtain a sensible estimate for the change in energy levels to 0(6) either by using 
perturbation theory or making a simple variational calculation. Since the effect is 
greatest for the ground state, we will analyse the displacement for the lowest energy 
level both ways - whose answers fortunately agree. 

To apply perturbation theory, take the unperturbed wave function in coordinate 
and momentum space: 

V7ra^ [pa^/n + ly ma 

and use both to work out the expectation value, 

AE = {b^Ry2m - ab/P) = ^ / {rtp{r)f dr - / p(j)\p) dp 

m Jo h"^ Jo 

= Zb'^a^/2m - IGaba/STrh ~ -16bn/37rm, (17) 

to first order in b. 

With the variational method adopt a trial wavefunction just like eq. (16), except 
that a is no longer identified with the Bohr radius. Then one gets (as a function of a) , 

, . fi^ 36^a^ a 16aba . . 

^^^^ = W + ^ ~ a " ^^^^ 
Minimising the energy, a cubic equation for the trial radius a is found: 

a IQab 



whose solution to order 6 is a ~ /ma + l%b%^ / 3m^ a'^ . Hence the optimal energy is 
obtained as 

E — —ma /2h^ - 16bh/3nm + ..., 
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thereby agreeing with (17). 

Taking the H-atom as the archetypical case we must make sure that atomic values 
are hardly touched by the reciprocity term ah/ P. The Rydberg energy of 13.6... eV 
~ 10"^*^ J is known to nine places of decimals and because /S.E ~ bh/m^ ~ 10~% J 
we shall require this to be < 10~^° J, just to be on the safe side. This sets a limit 
b < 10~^^ kg/s. This may seem miniscule but if we translate it into a universal angular 
frequency via Q = b/rrie, we obtain a value Q < 10^ Hz that is not so small. As we 
shall see presently, galactic considerations produce periods which are many orders of 
magnitude greater; therefore we can state with some confidence that the reciprocity 
modification with small enough b has essentially no effect on the accuracy of atomic 
calculations. 

4. Reciprocity and rotation of the galaxy 

We now treat the gravitational case, based on classical Newtonian dynamics and, 
presently, its reciprocity variant (3). It is a really important subject as the existence of 
possible dark matter owes a great deal to this study. In order to make any headway, 
observations of the motion/distribution of (visible) matter have to be taken into account 
[H H]. It has been established that the Milky Way contains a spinning supermassive 
black hole at centre with a mass of at least 4 x 10^ suns or about 10'^'^ kg. Visible mass is 
about 2 X 10^^ suns and on the basis of perceived rotation rates it is inferred that three 
or four times that mass is hidden in dark matter, as far out as we can see — assuming 
ordinary Newtonian gravity. The rotation speed of the Milky Way [5] seems to increase 
linearly with radius over a small distance, of about 0.3 kpc ~ 10^^ m, peaks, oscillates 
a bit and settles down to some 230 km/s out to 5 x lO^'' m from the black hole at the 
centre. We shall take this as a crude description of our galaxy, neglecting spiral arm 
structure and the effect of some barring of mass in the middle on the motion. We wish 
to investigate whether Born's reciprocity modification (3) has anything of consequence 
to say on the topic of dark matter and, in particular, if the observed tangential velocity 
profile is consistent with the visible matter distribution, without invoking dark matter. 

Since we are presuming that the rotation is mainly tangential, if we consider a 
typical star such as the sun, we may take — to a good approximation. So we just 
need to sum over all other masses M and their relative speeds v to obtain an average 
value for {a{l/r + b/p)) = {GmM{l/r + Q/v)). Here we have adopted a positive 6-sign, 
with b = rriQ > 0, as before, to fit in with the observations to come; the opposite sign 
gives nothing but grief. 

We cannot perform such averaging without some idea of the mass density 
distribution p(r). Unmodified Newtonian gravity indicates that the matter distribution 
is roughly spherical and p(r) ~ 1/r^ as far out as one can see, to ensure that the 
enclosed mass increases linearly with radius and reproduce a constant tangential star 
speed; hence dark matter. Of course the visible mass distribution fiatly contradicts this: 
it is mostly disk-shaped, with a concentration near the galaxy centre. This is what we 



Born Reciprocity & 1/r 



10 



shall model by taking a cylindrical Gaussian distribution, 

p{q, z) ^ /€25A<e~"'(^'+'^'^'Vf?7r', (20) 

to which we shall add a contribution from the central bulge including a black hole. 
(Here q is the distance from the central axis and z is the distance from the galactic 
plane.) The integral over the visible mass density produces a finite galactic mass M.\ 
the parameter k specifies the size of the galaxy disk while 5 is the ratio of disk diameter 
to disk thickness which varies from about 10 near the centre to 50 at the outer reaches 
of the disk; so we will take an average value 5^ ~ 1000 presently. 

Firstly we derive the potential energy in the plane of the disk at location {q, z = Q) 
for the cylindrical mass approximation above. 

V{g) = - ^^-^7"^ ^ / / / dcj>dz'dQ' , , , 

vr^ J J J ^ + (J^ — 2qq' cos (j) + z'^ 



TT 



J J dz'dg 



y/g'^ + g''^ + z' 



71 



GtiiMk f°° Q-UK g 

[21) 



TT 



where we have used the integral representation 2 ^ d^ = yrc/X. For large g note 

that V{g) ~ —GmAi/g as expected for a finite-sized source, while for small distances 
V{g) ~ 2GmM.hi^5g vanishes as ^ — > 0. To (21) we will shortly add a contribution from 
the central bulge, approximated by a mass cM. placed at the middle. 

The average over relative velocity f is a lot cruder and relies on data amassed over 
many years by astronomers. The galactic disk [3] exerts a strong attraction towards 
the plane on stars which stray from the disk and this results in a velocity dispersion 
Af^ ~ 30 km/s. As we are supposing that the majority of stars are turning at the same 
orbital speed w , the only relevant quantity is the angle between the azimuths of the two 



bodies and the velocity dispersion along the z-axis, |v' — v| = y (2f sin(0/2))2 + AAv^. 
Therefore the sum simplifies to an integration over azimuth: 



27r 1 r2TT 



27r|v'-v| 4:71V Jo ^^0/2+^ ttvVT+^ v 

where K is the elliptic integral of the first kind and e = Avz/v ~ 0.15 roughly [3]. (We 
use L as a parameter to fit the data in due course.) So without much compunction we 
shall simply take (fi/|v — v'|) ~ LQ/v as a fair approximation; after all we are only 
striving to get an idea of the size of the reciprocity constant Q here. 

Finally there is the matter of the halo contribution. Unlike mainstream ideas 
inferring dark matter, we rely on observations of visible matter (population II stars, 
white or brown dwarfs, star clusters, etc.) to put a bound of about ch — 10% on 
the amount of halo mass chM- Furthermore we shall suppose, like everyone else, that 
the halo velocity distribution is largely thermal [1]. In order to model these effects we 
will neglect the small oblateness of the halo and use a radial halo density distribution. 
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Pnir) ^ cj:/K/fA1/27r^r^(l + nj^r"^), multiplied by a normalized velocity distribution 
Ph{v) = (/3/7r)'^/^e~^"^. This fixes the mass of the halo to be and the halo velocity 

dispersion |4j to be {AvhY = 3/2/3, while the size of the halo is determined by 1/ un- 
it allows us to work out the contribution to the gravitational potential energy due to 
the halo: 

— / pH(r') 4:iTr''^dr' + / puir') Arrr' dr' 
r Jo Jr 



Vnir) = — Gm 



GmcHi^HM. 



TT 



2 1 
arctan(fi;j:/r) + ln(l H — ^ — -) 



(22) 



as well as the reciprocity halo contribution: 

= - / ph(v') A-nv'^dv' + / puW) A-nv'dv' = eiiU I3v)lv. (23) 

V V Jo Jv * 

Putting all this together, and including reciprocity terms, we obtain the total energy 
dependence on angular momentum i = rngv (conserved by axial symmetry) and distance 
from axis g for a test mass m: 



— + -mn^Q^-{L + c) --GmM - + ^ / 

'^2 I \ p Jtx jo 



du 



2mQ 
GmMcHi^H 



7r 



2 , X , , 1 

arctan(/t//r) + ln(l + 



^u{l + u){5^ + u)^ 



Gm'McH^Q ,,V^^, 
erf( ) 



Since we are presuming that the velocity is largely tangential, Pg = = dE/dg, which 
leads to the force equation: 

•2 n 



0= - 



mv 



mVt^Q - {L + c)GmM h 

VQ 

+ GmM Ht + 



+ 



\Q 

2GmMcH 



U 



n Jo ^ l + u){l + u/6^) 
GmMVtcH 



du 



arctan(fi;//f))-} 



2 ^^_^„._erf(vW 



TT 



.(24) 



We can simplify the look of this equation by rescaling to dimensionless variables. 
Let V = v/yjGMn, TZ = Kg, Lu = VL/\/GMh?,B = /3kGM. The tangential velocity 
profile then reduces to solving an equation for V as a function of TZ: 



uj(L + c) 
V' + ^ + chco 



erf(VgV) 
V 



2\l- e 

TT 



c ^ 27^ 
TZ \pk Jo 



ue 



du 



, 2ch fHH^l 

H — — arctan 



TXTI 



K 



.(25) 



'l + u){l + u/6^) 

The chosen sign for Q is highly significant in determining the behaviour of the velocity 
for small 71, viz. V ^ (L + c)TZuj /c, as required by the data. [Had we reversed the sign 
of or 6 we would not have been able to fit observations even remotely.] 
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Figure 4. Predicted tangential velocity curve with radius. 



A numerical solution of (23) is possible once a few measured values are input. 
The galaxy is 50000 light-years or more in radius, i.e. about 5 x 10^° m, so let us set 
the scale k ~ 0.7 x 10~^° m~^ or so. Taking the visible galactic mass to be roughly 
~ 2 X 10^^ kg, one estimates that k^GM ~ 5 x 10~^°/s^. Also the rotation speed 
outside the innermost part of the galaxy equals about 230 km/s and v grows to this value 
over about 0.3 kpc or 10^^ m, telling us that Q{L + c) / c ~ 2.5 x 10~^^ s~^. In attempting 
to fit the observed velocity profile, use values c ~ 0.05, as the proportion of galactic 
mass concentrated in the central bulge, and ch = 0.01 as the ratio of halo mass to disk 
mass. For simplicity take the halo radius to equal the disk radius or k = kh'-, although 
this is an underestimate it makes little difference to the numerical results because ch is 
quite small anyhow. A similar comment attaches to our chosen value of about 13 [1]. 
Finally set L ~ 0.85 and ~ 1.5 x 10~^^ Hz, implying uo ~ 0.2. None of these inputs 
is at all absurd. 

The profile equation (25) now determines V as a function of TZ and the numerical 
results are plotted in Figure 4. In trying to fit the data it is important to mention 
that there are two roots for the cubic in V. In the inner region we use the smaller 
root, corresponding to the linear relation v ~ rVL{L + c)/c, and beyond about r = 10^^ 
m we adopt the larger root. Although the Figure 4 graph is rather high below 2 kpc 
and shows a drop off to steadyish speed that is a bit faster than what is observed, the 
general shape of the curve is moderately satisfactory. Moreover the flatfish part of the 
curve has v ~ 200 km/s, which is the correct magnitude. However it must be admitted 
that our model and calculation have many rough edges and need a lot more refinement 
before they can be judged a success; it even conceivable that dressing the model may 
spoil it rather than enhance it. 
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5. Conclusions and criticisms 

The above fit has some good features and some bad ones; on the positive side it is able 
in principle to produce an orbital velocity curve that matches the data near the centre 
(linear rise) and at large distances (almost constant speed), without calling upon a 
dark matter component — Born's reciprocity principle has after all modified Newtonian 
dynamics in a radical way. But it must be confessed that the detailed fit above is 
far from convincing. This negative feature may be due to the crudeness of our visible 
mass and velocity distributions. It would be more realistic to include extra matter in 
the inner galactic bulge, introduce barring plus spiral arms and generally model the 
velocity distribution more accurately when working out the gravitational potential and 
its reciprocity counterpart. This is clearly fertile ground for future research. Meanwhile 
we can probably be content with our estimate offl ~ 1.5 x 10~^^ Hz for Born's hypothesis 
and not much more; it is truly a galactic scale and corresponds to a period of 1.3 x 10*^ 
years, which may be connected with the rotation rate of the galaxy. 

Serious concerns remain. The scheme destroys translational invariance because of 
the way that momentum and position are tied, so a violation of Galilean relativity is 
to be expected, never mind Lorentzian relativity. Choice of origin is another issue. We 
have picked the galactic centre as the obvious place; however one might fret about effects 
of harmonic acceleration (fi^r) which can get overwhelmingly large from outer regions 
of the universe, but is somewhat insignificant within a galaxy. Fortunately it appears 
that on cosmic scales the galaxies are uniformly distributed in every direction all around 
so one may presume that such attractions will cancel out overall. And because galaxies 
are receding away from us at Hubble rates the relative velocity term Q/v diminishes and 
balances out as one goes out. So it would seem, superficially at least, that one could 
choose any other galaxy and take the origin at its centre to reproduce its own galactic 
rotation without worrying unduly about other galaxies. 

Fundamental theoretical criticisms can nevertheless be levelled at (3). It is explicitly 
non-relativistic and should at least be made to conform with special relativity. To 
carry out that program sensibly one would need to study the quaplectic group [7] and 
augment the electromagnetic or gravitational field (photon/graviton exchange) by their 
reciprocity analogues or some other contributions. This signifies overhauling the whole 
of standard field theory and the task seems rather difficult, if not vague, at this stage. 
In the end it may turn out that Born reciprocity is unable to fit the data properly or 
mesh in with our familiar relativistic field theory concepts. Even if it works in limited 
fashion it would have to be seen as one of a panoply of Modified Newtonian Descriptions 
(MOND) of gravity over large distances. Taking a skeptical point of view, Born's idea 
will probably be found wanting, dark matter will be required and the problem of seeing 
it by non-gravitational methods will remain with us for a good while. But before the 
death knell is finally sounded on the subject of reciprocity, specialist galactic modellers 
need to investigate its ramifications comprehensively. 
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